Mon. Not. R. Astron. Soc. 000, 1-17 (2012) 



Printed 10 January 2013 (MN 1^1^ style file v2.2) 



A Bayesian approach to scaling relations for amplitudes of 
solar-like oscillations in Kepler stars 

E. Corsaro/'2'3* H.-E. Prohlich,^ A. Bonanno,^ D. Huber,^ T. R. Beddmg,*^'^ 
O. Benomar,^'^ J. De Ridder^, and D. Stello^'^ 

^Department of Physics and Astronomy, Astrophysics Section, University of Catania, Via S. Sofia 78, 1-95123 Catania, Italy 

'^I.N.A.F. - Astrophysical Observatory of Catania, Via S. Sofia 78, 1-95123 Catania, Italy 

^Instituut voor Sterrenkunde, K. U. Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium 

^Leibniz Institute for Astrophysics Potsdam (AIP), An der Sternwarte 16, 144^^ Potsdam, Germany 

^NASA-Ames Research Center, Moffett Field, CA 94035-0001, USA 

^Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, NSW 2006, Australia 
Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark 



Accepted . Received ; 



ABSTRACT 

We investigate different amplitude scaling relations adopted for the asteroseismology of 
stars that show solar-like oscillations. Amplitudes are among the most challenging astero- 
seismic quantities to handle because of the large uncertainties that arise in measuring the 
background level in the star's power spectrum. We present results computed by means of 
a Bayesian inference on a sample of 1640 stars observed with Kepler^ spanning from main 
sequence to red giant stars, for 12 models used for amplitude predictions and exploiting re- 
cently well-calibrated effective temperatures from SDSS photometry. We test the candidate 
amplitude scaling relations by means of a Bayesian model comparison. We find the model 
having a separate dependence upon the mass of the stars to be largely the most favored 
one. The differences among models and the differences seen in their free parameters from 
early to late phases of stellar evolution are also highlighted. 

Key words: methods: data analysis - methods: statistical - stars: oscillations - stars: 
solar-like - stars: late- type 



1 INTRODUCTION 

The study of solar-like oscillations (e. g. see Bedding k, Kjeld- 
sen 2003, 2008; Bedding 2011, for summaries and reviews) has 
experienced an enormous growth in the last five years thanks 
to the launch of the photometric space-based missions CoRoT 
(Baglin et al. 2006; Michel et al. 2008) and Kepler (Borucki 
et al. 2010; Koch et al. 2010). The latter in particular, is pro- 
viding a very large amount of high quality light curves, with a 
very high duty cycle (see Gilliland et al. 2010b; Brown et al. 
2011; Garcia et al. 2011, for a general introduction to the Ke- 
pler asteroseismic program, for the presentation of the Kepler 
Input Catalog, and for a description of the preparation of Ke- 
pler light curves, respectively). These will become longer with 
the upcoming era of the Kepler extended mission (Still 2012). 

* E-mail: eco@oact.inaf.it; 



These asteroseismic studies have recently led to the birth 
of the ensemble asteroseismology (Chaplin et al. 2011a), which 
is showing great potential for a thorough understanding of stel- 
lar evolution theory. The success of ensemble asteroseismology 
relies mainly on adopting scaling relations: generally simple em- 
pirical laws that allow for the derivation of fundamental stellar 
properties for stars different from the Sun by scaling their as- 
teroseismic quantities from the solar values. 

Among the most challenging asteroseismic quantities to 
measure and model one can certainly mention the oscillation 
amplitude. This is due to both the difficulty in estimating the 
background level in the power spectrum and the rather compli- 
cated physics involved in the driving and damping mechanisms 
of the modes (e. g. see Kjeldsen et al. 2005, 2008). Different scal- 
ing relations aimed at predicting amplitudes by scaling from 
the Sun's values have been derived and discussed by several 
authors, both theoretically (Kjeldsen & Bedding 1995; Houdek 
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et al. 1999; Houdek & Gough 2002; Houdek 2006; Samadi et al. 
2007; Belkacem et al. 2011; Kjeldsen & Bedding 2011; Samadi 
et al. 2012) and observationally (Stello et al. 2011; Huber et al. 
2011b; Mosser et al. 2012). A variety of amplitude scaling rela- 
tions has been used extensively in literature in both ensemble 
studies (e.g. see Stello et al. 2010; Chaplin et al. 2011b; Verner 
et al. 2011; Huber et al. 2011b; Mosser et al. 2012; Belkacem 
2012) and detailed analyses of single stars from main sequence 
to the subgiant phase of evolution (e. g. Bonanno et al. 2008; 
Samadi et al. 2010; Huber et al. 2011a; Mathur et al. 2011; 
Campante et al. 2011; Corsaro et al. 2012a). 

The underlying physical meaning of these various ampli- 
tude scaling relations is still not properly understood (e. g. see 
the discussions by Samadi et al. 2007; Verner et al. 2011; Hu- 
ber et al. 2011b; Samadi et al. 2012; Belkacem 2012). Testing 
them with observational data is vital for assessing the com- 
peting relations and for improving our understanding of stellar 
oscillations, i.e. the driving and damping mechanisms that pro- 
duce the observed amplitudes (see also Chaplin et al. 2011b). 
In this context, Bayesian methods can be of great use (see e. g. 
Trotta 2008) because they allow us to measure physical quanti- 
ties of interest in a rigorous manner. Moreover, Bayesian statis- 
tics provides an efficient solution to the problem of model com- 
parison, which is the most important feature of the Bayesian 
approach (see also Benomar et al. 2009; Handberg & Campante 
2011; Gruberbauer et al. 2012). 

In the present paper we analyze amplitude measurements 
of a sample of 1640 stars observed with Kepler, together with 
temperature estimates derived from SDSS photometry, which 
we introduce in Section 2. In Section 3 we discuss the differ- 
ent scaling relations for predicting the oscillation amplitude 
per radial mode. The results obtained from a Bayesian param- 
eter estimation for the different scaling relations are shown in 
Section 4 and the model comparison is presented in Section 5. 
Lastly, discussion and conclusions about the results of our anal- 
ysis are drawn in Section 6 and 7, respectively. 



2 OBSERVATIONS AND DATA 

We use amplitude measurements and their uncertainties, ob- 
tained by Huber et al. (2011b) for a sample of 1673 stars span- 
ning from main sequence (MS) to red giant stars (RGs) ob- 
served with Kepler in short cadence (SC; mostly MS stars but 
also some subgiants and low luminosity RGs) and long cadence 
(LC; ah RGs) modes (Gilliland et al. 2010a; Jenkins et al. 2010, 
respectively). Most of the 542 stars observed in SC have pho- 
tometry for one month, while the 1131 stars observed in LC 
have light curves spanning from Kepler^s observing quarters 
to 6. All amplitudes were derived according to the method 
described by Kjeldsen et al. (2005, 2008), which provides ampli- 
tudes per radial mode (see Huber et al. 2009, for more details). 
Values of the frequency of maximum power, i^max, the large fre- 
quency separation, Ai/, and their uncertainties for all the stars 
were also taken from Huber et al. (2011b), who used the SYD 
pipeline (Huber et al. 2009). 

It is important to have accurate temperature estimates for 



the stars of our sample. Unfortunately, those provided by the 
KIC (Brown et al. 2011) are known to suffer from significant 
systematic effects (see Pinsonneault et al. 2012a for a detailed 
discussion of the problem). We used revised Tgff derived by 
Pinsonneault et al. (2012a) for a total of 161977 KIC stars 
from Sloan Digital Sky Survey griz filters, which were corrected 
using temperature estimates from infrared flux method (IRFM) 
{J — Kg) color index for hot stars (e.g. see Casagrande et al. 
2010). The revised effective temperatures are available in the 
online catalog (Pinsonneault et al. 2012b). 

By cross-matching the stars of our sample with the tem- 
perature estimates provided by Pinsonneault et al. (2012a), we 
arrived at a final sample of 1640 stars with an accurate Teff 
(1111 observed in LC and 529 in SC), which will be used for our 
investigation. Total uncertainties on temperature, as derived by 
Pinsonneault et al. (2012a), include both random and system- 
atic contributions. The amplitudes of the final sample are plot- 
ted against z^max and A in in Figure 1 (top and middle panels, 
respectively) . The bottom panel shows an asteroseismic HR di- 
agram of our sample of stars (amplitudes against Teff from Pin- 
sonneault et al. 2012a), similar to the one introduced by Stello 
et al. (2010). la error bars are overlaid on both quantities for 
each panel. The average relative uncertainty in amplitude is 
{(Ja/A) = 9.2% for the entire sample, and (crAM)sc = 11.7% 
and ((Ja/^)lc = 8.1 %, for SC and LC targets, respectively. 



3 AMPLITUDE SCALING RELATIONS 

Several scaling relations for oscillation amplitudes have been 
proposed so far. We will briefly introduce them in the following 
(see Huber et al. 2011b for further discussion). 



3.1 The L/M scaling relation 

The first scaling relation for amplitudes was introduced by 
Kjeldsen &; Bedding (1995) for radial velocities, based on the- 
oretical models by Christensen-Dalsgaard &; Frandsen (1983). 
It is given by 



(1) 



where Vosc represents the prediction for the amplitude in radial 
velocity, L is the luminosity and M the mass of the star, and 
s = 1. Kjeldsen & Bedding (1995) also showed that the corre- 
sponding photometric amplitude Ax , observed at a wavelength 
A, is related to Vosc by 



Ax oc 



(2) 



For adiabatic oscillations, the exponent r is 1.5 (see also Houdek 

2006) , but Kjeldsen & Bedding (1995) found the observed value 
for classical p-mode pulsators to be 2.0 (see also Samadi et al. 

2007) . By combining the two equations, one obtains 



Ax oc 



(3) 
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Figure 1. Oscillation amplitudes for 1640 stars observed with Ke- 
pler in SC (orange squares) and LC (blue circles) modes and plotted 
against the frequency of maximum power i^max (top) and the large 
frequency separation Au (middle) of the stars in a log-log scale. Am- 
plitudes against the effective temperature Teff are shown in the bot- 
tom panel, representing an asteroseismic HR diagram for our sample 
of stars. Icr error bars are shown on both quantities for all the plots. 
The Sun is shown with its usual symbol (0). 



We are interested in a sample of stars observed with Kepler, 
whose bandpass has a central wavelength A = 650 nm. By scal- 
ing Eq. (3) to our Sun, we have 



\M/Mq 



(4) 



where Aeso,© = 3.98 ppm is the Sun's photometric amplitude 
observed at the Kepler wavelength (e.g. see Stello et al. 2011) 
and Teff,0 = 5777 K. The exponent s has been examined both 
theoretically (e.g. see Houdek et al. 1999; Houdek & Gough 
2002; Samadi et al. 2007) and observationally (Gilliland 2008; 
Dziembowski & Soszynski 2010; Stello et al. 2010; Verner et 
al. 2011; Baudin et al. 2011), and found to be 0.7 < s < 1.5. 
The exponent r has usually been chosen to be 2.0 (e.g. see 
the discussion by Stello et al. 2011), which implies that solar- 
like oscillations are not fully adiabatic. This was also shown 
in the case of CoRoT red giant stars by Samadi et al. (2012). 
Nonetheless, some authors (e. g. see Michel et al. 2008; Mosser 
et al. 2010) have chosen to adopt r = 1.5 (see also the discussion 
by Kjeldsen & Bedding 2011). 

We can also use the results introduced by Brown et al. 
(1991), which suggest that the frequency of maximum power 
scales with the cut-off frequency of the star. Hence, z^max oc 
g / V^eff, and by considering that L/M scales as T^^/g, one 
obtains 



(5) 



TT 

M Vm 

Thus, Eq. (4) can be rewritten as 
Ax 



teff 



3.5s- 



(6) 



^A,0 V^'max,© / V^eff,0 

where z^max,© = 3100 /iHz. The functional form of the ampli- 
tude scaling relation given by Eq. (6) has the advantage of 
simplifying the inference problem presented in Section 4 with 
respect to the one of Eq (4). This is because Eq. (6) is based on 
a set of independent observables, namely z^max and Teff, which 
represent the input data used for this work. An analogous argu- 
ment has been applied to the other scaling relations described 
in the following sections. 

Thus, Eq. (6) represents the first model to be investigated, 
which we will refer to as Mi. For this model, both the ex- 
ponents s and r are set to be free parameters. In parallel, an 
extended version of Eq. (6) given by 

^A,© V^max,©/ V^eff,©/ 

is also considered, where the factor f3 allows the model not to 
necessarily pass through the solar point, as we will discuss in 
more detail in Section 4.2. Eq. (7) is treated as a separate model 
and will be denoted as Mi^p. Clearly, Mi,p depends upon the 
additional free parameter, /3. 



3.2 The bolometric amplitude 

The stars considered span over a wide range of temperatures, 
from about 4000 K to more than 7000 K. Hence, following the 
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discussion by Huber et al. (2011b), a more valuable expres- 
sion for the photometric amplitude could be represented by 
the bolometric amplitude Aboi, which is related to Ax by (see 
Kjeldsen & Bedding 1995): 



Abol (X XAx Teff oc 



By using Eq (1) we thus have 



i(i) 



M 



which by scaling to the Sun and adopting Eq. (5) yields 



A 



(1) 



^bol,0 



^max,0 



3.5s-r+l 



(8) 



(9) 



(10) 



where Aboi,© = 3.6 ppm represents the Sun's bolometric am- 
plitude, determined by Michel et al. (2009) and also adopted 
by Huber et al. (2011b). Eq. (10) is the second model, M2, 
to be investigated in Section 4, with the exponents s and r 
set again to be the free parameters. We also consider the new 
model M.2,/3, which again includes the proportionality term /3 
playing the same role as in Eq. (7). 



3.3 The Kepler bandpass-corrected amplitude 

Ballot et al. (2011) have recently established a bolometric cor- 
rection for amplitude of radial modes observed with Kepler, 
which translates into a correction factor for effective tempera- 
tures falling within the range 4000-7500 K. Again following the 
approach by Huber et al. (2011b), we consider a revised version 
of Eq. (9), which reads 



l(2) 



^eff ^CK(Teff) 



where 



Ck {Tei 



1 5934/ 



(11) 



(12) 



is the bolometric correction expressed as a power law of the 
effective temperature. By scaling once again to the Sun and 
applying Eq. (5), we obtain 



l(2) 



Ab, 



ol,0 



^max,0 



teff,0 



3.5s-r+0.2 



(13) 



which we will refer to as M.3- As for the other models, we intro- 
duce the model A^3,/3 with the proportionality term /3 included. 



3.4 The mass-dependent scaHng relation 

A mass dependence of the oscillation amplitudes was suggested 
for the first time by Huber et al. (2010), and later on studied in 
detail by Stello et al. (2011) for cluster RGs with the introduc- 
tion of a new scaling relation. It was also tested afterwards by 
Huber et al. (2011b) for a wider sample of field stars. According 



to Huber et al. (2011b), an obvious way to modify Eq. (11) is 
given by 



^bol ^ 



(14) 



where now the mass varies with the independent exponent t. 
For this case, the dependence upon the quantities z^max and Tgff 
becomes slightly more complicated because the simple propor- 
tionality expressed by Eq. (5) can no longer be adopted. There- 
fore, the first step to derive the functional form based on our 
set of observables (z^max, Az/, Tgff), is to consider the scaling re- 
lations for the large frequency separation Av (e. g. see Bedding 
et al. 2007) 



R 

Rq 



(15) 



with Az/0 = 135 /iHz, and for the frequency of maximum power 



Z^max,0 



M 

Mo 



R 

Rq 



(16) 



both expressed in terms of the fundamental stellar properties 
M, R, and Tgff. It is however worth mentioning that these scal- 
ing relations are empirical approximations whose validity and 
limitations are not yet fully understood and currently under 
debate (e.g. see Miglio et al. 2012, and references therein). By 
combining Eq. (15) and (16), one can derive an expression for 
the seismic radius of a star (e.g. see Chaplin et al. 2011a), 
namely 



R 

Rq 



Z^max,0 



v 



Avq) \TeH,Q 



(17) 



As a second step, we express /M* in terms of R, z^max, and 
Teff and scale to the Sun's values, yielding 



L 



M 
Mo 



R 

Rq 

Teff 
^eff,0 



Z^max,0 
4s-0.5t 



(18) 



Finally, by combining Eqs. (12), (14), (17) and (18) we arrive 
at 



l(3) 



^bol,( 



Z^max,0 



(—Y~ 



-1.5t-r+0.2 



(19) 



This represents the model for the mass-dependent scaling rela- 
tion for amplitudes, hereafter denoted as In this case, we 
have the three free parameters s, r, and t and the set of ob- 
servables now includes also our measurements of Av. The cor- 
responding model A44,/3 has the largest number of free param- 
eters among those investigated in this work. Note that models 
M4 and A44,/3 reduce to models M.3 and A^3,/3, respectively, 
for t = s. 
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3.5 The lifetime-dependent scaling relation 

Kjeldsen & Bedding (2011) have recently provided physical ar- 
guments to propose a new scaling relation for predicting the 
amplitudes of solar-like oscillations observed in radial veloci- 
ties. Their relation arises by postulating that the amplitudes 
depend on both the stochastic excitation (given by the gran- 
ulation power, see Kjeldsen & Bedding 2011, for details) and 
the damping rate (given by the mode lifetime). It reads 



(20) 



where Tqsc is the average mode lifetime of radial modes. By 
means of Eq. (8), and with the bolometric correction introduced 
in Section 3.3, the corresponding relation for the bolometric 
amplitude is given by (see also Huber et al. 2011b) 



' osc 



(21) 



In order to obtain the expression for the model we use 

similar arguments to those adopted in Section 3.4, arriving at 



l(4) 
^bol 



ibol,0 



^max,0 



'T'osc,© 



(22) 



where Tosc,© = 2.88 d, as adopted by Kjeldsen & Bedding 
(2011). For our computations we assume that the mode life- 
time is a function of the effective temperature of the star alone 
(e.g. see Chaplin et al. 2009; Baudin et al. 2011; Appourchaux 
et al. 2012; Belkacem et al. 2012; Corsaro et al. 2012b). We 
used the empirical law found by Corsaro et al. (2012b), which 
relates the mode linewidths F of the radial modes {£ — 0) to the 
effective temperature of the star within the range 4000-7000 K. 
In particular, they found that 



F = Fo exp 



To 



(23) 



where Fq = 1.39 ± 0.10 /iHz and To = 601 ± 3K. This relation 
was calibrated using Kepler RGs in the open clusters NGC 6791 
and NGC 6819, and a sample of MS and subgiant Kepler field 
stars. Given r = (ttF)"^, we obtain 



Tosc = To exp 



?eff,0 ~ ?eff 



(24) 



with To = 2.65 =b 0.19 d. The mode lifetimes were computed 
for all the stars of our sample by means of Eq. (24), together 
with their corresponding uncertainties from the error propaga- 
tion. As for the other scaling relations, we also introduce model 



3.6 A new scaling relation 

Following similar arguments to those adopted by Stello et al. 
(2011) for introducing a new scaling relation for amplitudes of 
cluster RGs, and the discussion by Huber et al. (2011b) about 



the mass dependence of the amplitudes, we introduce a slightly 
modified version of the amplitude relation given by Eq. (21), 
where we set the mass to vary with an independent exponent 
t, thus yielding 



A 



(5) 



Tli'+^CK{T^^)' ^^^^ 
By adopting again Eq. (18) and rearranging, we finally obtain 



l(5) 



^bol,0 



^max,0 
Teff 



f— V 



Tosc,0 



4.55-r-1.5t 



(26) 



^Teff,0 ^ 

hereafter marked as model Mq. A.s done for the other scaling 
relations, the model Mg,^ is also included in our inference. 
Clearly, models Me and M6,/3 reduce to models M5 and A^5,/3, 
respectively, for t = 1.5. 



4 BAYESIAN INFERENCE 

We now use Bayesian inference for the free parameters of the 
models described above. The Bayes' theorem tells us that 

p{A\i,MMi\M) 



P{$\A,M) 



(27) 



p{A I M) 

where ^ = ^i,^2,...,^fc is the vector of the k free parameters 
that formalize the hypotheses of the model considered for 
the inference, and A is the set of amplitude measurements. The 
term p{A \ ^, Ai) is now identified with the likelihood C{^) of 
the parameters given the measured oscillation amplitudes: 



piA\^,M) = /:i^\A,M). 



(28) 



Thus, the left-hand side of Eq. (27) is the posterior probabil- 
ity density function (PDF), while the right-hand side is the 
product of the likelihood function £(^), which represents our 
manner of comparing the data to the predictions by the model, 
and the prior PDF 7r(^ | M), which represents our knowledge of 
the inferred parameters before any information from the data is 
available. The term p{A \ J\4) is a normalization factor, known 
as the Bayesian evidence, which we do not consider for the in- 
ference problem because it is a constant for a model alone. As 
we will argue in Section 5, the Bayesian evidence is essential 
for solving the problem of model comparison. 

For our inference problem, we adopt the common Gaussian 
likelihood function, which presumes that the residuals arising 
from the difference between observed and predicted logarithms 
of the amplitudes are Gaussian distributed, i. e. the amplitudes 
themselves are presumed to be log-normal distributed (see also 
Appourchaux et al. 1998). Therefore, we have 



'^(^) = 11 ^''P ~o[~^ 

V27ro-i [ 2 V 



(29) 



where N is the total number of data points (the number of 
stars, in our case), while 

Ai(0 = ln^f^-ln4''($) (30) 
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is the difference between the observed logarithmic amphtude 
for the i-th star and the predicted one (which depends on the 
adopted model, i.e. on the parameters vector ^). The term 
ai appearing in the leading exponential term of Eq. (29) is 
the total uncertainty in the predicted logarithmic amplitude, 
namely the relative uncertainty of the amplitude enlarged by 
the relative errors of the independent variables i^max, Ai/, and 
Teff . This means that we are not considering error- free variables 
in our models (see e. g. D'Agostini 2005; Andreon & Hurn 2012, 
for more details). For simplifying the computations a modified 
version of the likelihood function, known as the log-likelihood, 
is preferred. The log-likelihood function is defined as A(^) = 
ln£(^), which yields 



MO 



where 



1 ^ 



In V27T(7i . 



(31) 



(32) 



The choice of reliable priors is important in the Bayesian ap- 
proach. For our purpose, uniform priors represent a sensible 
choice for most of the free parameters. This means one has no 
assumptions about the inferred parameters before any knowl- 
edge coming from the data, with equal weight being given to 
all values of each of the parameters considered. In particular, 
we use standard uniform priors for the exponents s, r and t of 
the models described above, letting the parameters vary within 
a limited range in order to make the priors proper, i. e. normal- 
izable to unity. For the proportionality term /3 introduced in 
Eq. (7), we adopt the Jeffreys' prior oc /3~'^ (Kass & Wasser- 
man 1996), a class of uninformative prior that results in a uni- 
form prior for the natural logarithm of the parameter. In this 
manner, the parameter of interest is represented by the offset 
\n/3 (see below), whose prior is uniform distributed and also 
limited in range. Hence, uniform priors are included in the in- 
ference problem as a simple constant term depending on the 
intervals adopted for the inferred parameters 



4 



max 
j 



(33) 



with a = ^2 = ^3 = t and ^4 = In^, and 
the minimum and maximum values defining the interval of the 
j-th parameter. The intervals that we adopt are listed in Ta- 
ble 1. These ranges are used for both the Bayesian parameter 
estimation and the model comparison. 

A note of caution concerns the treatment of the uncer- 
tainties. In fact, by using the natural logarithm of the equa- 
tions that describe the models (see also Section 4.2), we ensure 
that we are not favoring for example frequencies upon periods, 
amplitudes upon power, temperatures upon surface brightness, 
etc., which has the advantage of making the error propagation 
law fully correct. Thus, the corresponding uncertainties to be 
considered in Eq. (29) will be the relative uncertainties. 

The inference problem for a given parameter, e.g. ^1, is 



Table 1. Maximum and minimum values of the free parameters, 
adopted for all the models and samples. 





tmin tmax 




S 


[0.2,1.2] 


r 


[-6.5,11.0] 


t 


[1.0,2.0] 


ln/3 


[-1.0,1.0] 



then performed by integrating (marginalizing) the posterior 
distribution function \ d,M) over the remaining k — 1 pa- 
rameters ^2, ^3, • • • , ^fc- We obtain the corresponding marginal 
PDF of the parameter ^1 



pi^i\A,M) 



I A,M)d£,2,d£.3,...,dS,k, 



(34) 



whose statistical moments and credible intervals (i. e. Bayesian 
confidence intervals) are the quantities of interest for our work. 
Since the dimensionality of our problem is not higher than 
= 4, all the integrations can be computed by direct numerical 
summation of the posterior distribution over the remaining pa- 
rameters (only in case the observables, like Tgff , were error-free 
all integrations could be computed analytically). 

The results presented in the coming sections are derived 
in three cases: for the entire sample of stars, and for SC (domi- 
nated by MS stars) and LC (RGs) targets separately. Analyzing 
the two subsets separately allows us to test whether the fitted 
parameters of the scaling relations are sensitive to the evolu- 
tionary stage of the stars (see also Huber et al. 2011b). 

The mean values (or expectation values) of the free param- 
eters of the models, together with their corresponding 68.3 % 
Bayesian credible intervals, are listed in Table 2 for the case of 
the entire sample, and in Tables 3 and 4 for the subsets of SC 
and LC targets. We also computed a weighted Gaussian rms, 
c^^?ns5 of the residuals Af(^), where we adopted the weights 
Wi = d-^^, ai being the total uncertainty used in Eq. (29). 
The maximum of the log- likelihood function, Amax, and crj^^ns? 
used as an estimate of the fit quality, are also listed in the 
same tables. In addition, we derived correlation coefficients for 
each pair of free parameters by means of principal component 
analysis using Singular Value Decomposition (SVD) from the 
posterior PDFs. The results are shown in Tables 5, 6, 7 for 
the entire sample, and SC and LC targets, respectively, with 
— 1 meaning total anti-correlation and 1 total correlation. The 
effects of the correlations will be discussed in Section 6. 



4.1 Statistically independent models 

Before getting to a description on how to correctly include the 
models in the inference problem, it is useful to highlight that 
the models A4i,/3, A42,/3, and Ads,^ on the one hand, and the 
models A^2, and M.3 the other hand, are statistically identical 
to one another (but not identical in the general sense, since 
their underlying physical assumptions are different). In case 
the intervals of their free parameters are the same for all the 
models when performing the Bayesian parameter estimation (as 
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it is for the analysis presented in this work) this imphes that a 
statistical inference for these models would lead to identical (or 
directly related) values of these free parameters. In particular, 
according to Eqs. (7), (10), and (13) we have that 

s(A^l,/3) = s(A^2,/3) = s(A^3,/3) 

s{M2) = siMs) 
for the exponent s, 

r{Mi,i3) = r{M2,f3) - 1 = r{M3,p) - 0.2 
r{M2) = riMs) - O.S , 
for the exponent r, and 

^bol,0 



(35) 



(36) 



ln/3(A4i,^) = ln/3(A42,/3) + ln 
= ln/3(A^3,/3)+ln 



^A,0 
^bol,0 
^A,0 



(37) 



for the offset ln/3, where the term In (Aboi,0/^A,0) arises from 
the difference in considering the amplitudes to be either ob- 
served at A = 650 nm (A^i,/?) or bolometric {M.2,/3 and A^s,^)- 
As a consequence, from now on the entire analysis for the mod- 
els A^i,/3, Al2,/3, and Ms,^ on one side, and for the models M2 
and M3 on the other side, will be reduced to that of the two 
models M.i,p and AI2, respectively. The reader can derive the 
corresponding parameters for the other dependent models using 
Eqs. (35), (36), and (37). 



4.2 Models Mi and Mi,/3 

The first model to be investigated is given by Eq. (6). As argued 
above, we need to consider the natural logarithm in order to 
treat the observables independently of the function adopted 
(see the discussion in Section 4). Hence the model reads 



In 



= — s In 



^max,0 



+ (3.5s - r) In 



(38) 



At this stage we briefly describe how the uncertainties have 
been included in our analysis. The new uncertainties on the 
scaled amplitude to be considered in Eq. (29) are clearly given 
by aAi/Ai, hereafter for simplicity, where Ai is the observed 
amplitude for the i-th star and cta^ its corresponding uncer- 
tainty as derived by Huber et al. (2011b). However, Eq. (38) 
suggests that the uncertainty on amplitude is not the only one 
affecting the predicted amplitude Ax. In fact, uncertainties on 
both z^max (derived by Huber et al. 2011b) and Tgff (from Pin- 
sonneault et al. 2012a) have to be included in our computa- 
tions. The total uncertainty to be used in Eq. (29) is given by 
the Gaussian error propagation law, which gives 



(Ji (s,r) = aA^-\- s (7^^ 



+ {?^.^s-Tfal^^^. (39) 



where we defined cft^ 



and CTTeff^i = 



^^eff i/^eff,i, similarly to what was done for the amplitudes. 



Eq. (39) only holds in case of uncorrelated uncertainties and 
linear relations (see also D'Agostini 2005; Andreon & Hurn 
2012). Intuitively, Eq. (39) is the quadratic sum of the relative 
uncertainties over the physical quantities considered, according 
to Eq. (38). A variation of Eq. (38) is represented by the model 
A4i,/3 given by Eq. (7), whose natural logarithm reads 



In 



Ax 



A 



A,0 



s In 



^max,0 



-h (3.5s - r) In 



^eff,0 



(40) 



+ ln/ 



which differs from Eq. (38) by the additional term In /3. As al- 
ready argued before, the offset ln/3 allows the model not to 
necessarily pass through the solar point (A©, z/max,0, Teff,©). 
Its introduction in the inference is of importance if one wants 
to assess whether or not the Sun is a good reference star for 
the sample considered. This choice is also motivated by the 
fact that the Sun is falling at the edge of the sample of stars 
when plotting amplitudes against z^max and Azy (Figure 1, top 
and middle panels). This peculiar position is also evident from 
our asteroseismic HR diagram (Figure 1, bottom panel), and 
is caused by the lack of solar twins in our sample of stars (see 
the discussion by Chaplin et al. 2011a). In fact, in the case of 
In /3 7^ 0, by replacing the solar values in Eq. (40) (or alterna- 
tively Eq. (7)), the predicted amplitude for the Sun would be 
resized by a factor (3. This means that the best reference ampli- 
tude for scaling the amplitudes of our sample of stars would be 
represented by I3Ax,q. According to Eq. (40), the total uncer- 
tainty to be considered in building the likelihood function for 
the model Mi,j3 is given again by Eq. (39) because the offset 
does not play any role in the total contribute of the uncertain- 
ties. 

A representative sample of the resulting marginal PDFs 
is plotted in Figure 2 for the case of the entire sample, where 
68.3% Bayesian credible regions (shaded bands) and expec- 
tation values (dashed lines) are also marked. The comparison 
between the predicted and the observed amplitudes is shown in 
Figure 3 for the three cases considered (top panels for model 
All, bottom panels for model A4i,/3), together with a plot of 
the residuals arising from the difference between the models 
and the observations. 



4.3 Model M2 

The model M2 given by Eq. (10) deserves a similar description 
to that presented in Section (4.2) for model Ali, where the 
natural logarithm is now given by 



4(1) 

In' 



A, 



bol,0 



s In 



^max,0 



-h (3.5s - r -h 1) In 



(41) 



^7'eff,0 ^ 

with a total uncertainty for the i-th star expressed as 

a- (s, r) = a\. + s^al^^^ . + (3.5s -r+ if'^l^,,^^ (42) 
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Table 2. Expectation values of the inferred parameters for all the models described in Section 3 in the case of the entire sample (both LC and 
SC targets), having N = 1640 stars. 68.3% Bayesian credible intervals are added. The maximum value for the log- likelihood function Amax 
and a weighted Gaussian rms of the residuals, crj^ns? ^^^^ reported. 



Model 


s 


r 


t 


In^ 


A-max 


w 


Ml 


0.680 ± 0.002 


4.31 ± 0.04 






-3533.1 


0.23 


Mi^p 


0.524 di 0.004 


5.51 di 0.05 




0.400 ± 0.010 


-2491.3 


0.24 


M2 


0.722 di 0.002 


4.96 =b 0.04 






-4159.9 


0.24 


M4 


0.822 dz 0.003 


3.93 lb 0.06 


1.58 ±0.02 




-948.5 


0.23 


M4,/3 


0.643 di 0.005 


4.46 ± 0.06 


1.36 ±0.02 


0.528 ±0.012 


163.4 


0.18 


M5 




-5.71 ±0.03 






-925.5 


0.31 


M5,/3 




-5.07 di 0.05 




-0.122 ±0.007 


-786.9 


0.29 


Me 




-5.04 ±0.05 


1.80 ±0.02 




-727.1 


0.30 


M6,/3 




-4.97 ±0.06 


1.75 ± 0.02 


-0.035 ± 0.012 


-722.8 


0.29 



Table 3. Same description as for Table 2 but in the case of SC targets, having N = 529 stars. 



Model s r t h\(3 Amax cr\ 



Ml 


0.775 ± 0.003 


3.23 ± 0.05 






-321.4 


0.20 


Mi^p 


0.624 ±0.010 


3.68 ± 0.06 




0.241 ±0.016 


-199.3 


0.19 


M2 


0.838 ± 0.003 


4.05 ±0.05 






-443.8 


0.22 


M4 


U.y84_Q Q^Q 


2.79 ± 0.09 


l-66_o.o5 




-94.9 


0.23 


M4,/3 


0.748 ±0.015 


3.47 ±0.09 


1.27 ±0.04 


0.321 ± 0.020 


18.0 


0.20 


Mb 




-2.78 ±0.09 






-83.1 


0.24 


M5,/3 




-2.75 ± 0.09 




0.020 ± 0.008 


-80.2 


0.24 


Me 




9 on+0-10 
-^.»u_Q Q9 


1.56 ±0.02 




-79.6 


0.24 


Me, 13 




-2.75 ±0.10 


1.72 ±0.04 




-59.4 


0.24 



Table 4. Same description as for Table 2 but in the case of LC targets, having N = 1111 stars. 



Model 


s 


r 


t 


In^ 


Amax 


w 


Ml 

Mi^p 

M2 

Ma 
M4,p 


464+°-°°^ 
0.548 ± 0.009 

U.^yi_Q QQg 

U.DDD_Q QQ5 

0.602 ± 0.008 


y-O<J_0.l6 

9.671^1^6 

10.511^1^6 
6.99 ±0.12 
5.87 ±0.14 


1.28 ±0.02 
1.31 ±0.02 


-0.35 ±0.03 
0.45 ± 0.03 


-799.2 
-737.4 
-769.0 
207.4 
301.0 


0.28 
0.27 
0.28 
0.18 
0.16 


Ms 




-6.08 ± 0.04 






-357.9 


0.24 


M5,p 
Me 




-4.38 ±0.16 
-5.94 ±0.08 


1.55 ±0.03 


-0 27+°-°^ 
-0.03 


-309.2 
-356.0 


0.24 
0.24 


Me,^ 




-4.39 ±0.16 


1 45+0-02 


-0.30 ±0.03 


-307.1 


0.24 



Table 5. Correlation coefficients for pairs of free parameters for each model in the case of the entire sample. 



Model 


s vs r 


s vs 6 


s vs t 


r vs 6 


r vs t 


6 vs t 


Ml 


-0.90 












Mi^l3 


-0.85 


-0.92 




0.66 






M2 


-0.94 












M4 


-0.74 




0.59 




0.01 




M4,(3 


-0.74 


-0.76 


0.44 


0.31 


-0.25 


0.04 


Mb 














Mb,^ 








-0.74 






Me 










0.74 




Me,/3 








-0.41 


0.18 


0.71 
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Table 6. Same description as for Table 5 but in the case of SC targets. 



Model 


s vs r 


s vs 6 


s vs t 


r vs 6 


r vs t 


6 vs t 


Ml 


0.22 














-0.40 


-0.94 




0.49 






M2 


0.20 












M4 


-0.29 




0.88 




-0.27 




M4,p 


-0.48 


-0.85 


0.71 


0.44 


-0.37 


-0.31 


M5 






















0.12 






Me 










-0.07 




M6,/3 








0.03 


-0.01 


0.76 



Table 7. Same description as for Table 5 but in the case of LC targets. 



Model 


s vs r 


s vs 6 


s vs t 


r vs 6 


r vs t 


6 vs t 


Ml 


-0.98 












Ml, 13 


-0.49 


-0.71 




-0.25 






M2 


-0.98 












Ma 


-0.92 




0.50 




-0.19 




M4,^ 


-0.31 


-0.66 


0.30 


-0.46 


-0.30 


0.18 


Mb 














M5,/3 








-0.98 






Me 










0.88 




Mq^P 








-0.88 


0.04 


0.37 



to be included in Eq. (29). The resulting models are shown in 
Figure 4, with similar descriptions as those adopted for Fig- 
ure 3. 



4.4 Models Ma and Ma,^ 

The models Ma and Ma,i3 (see Section 3.4) are clearly the 
most complex ones among those investigated in this work be- 
cause the largest number of free parameters is involved, and 
measurements of Azy are also needed. We note that, although 
a tight correlation between i^max and Ai/ has been found in 
previous studies (e. g. see Stello et al. 2009a), we choose not to 
adopt the z^max-Azy relation to express model Ma in terms of 
z^max only (or alternatively Az/) because additional uncertain- 
ties arising from the scatter around this relation would affect 
the results of our inference. This is also motivated by recent 
results by Huber et al. (2011b) who found that the 
relation changes as a function of Teff between dwarf and giant 
stars. 

Therefore, by considering the natural logarithm of Eq. (19) 
one obtains 



In 



l(3) 



^bol,G 



(2s - 3t) In 



Z^max,0 



+ (5s- 1.5t-r + 0.2)ln 



^eff,0 



(43) 



for model Ma^ and with the additional term In^ for Ma.ij 



According to Eq. (43), the total uncertainty to be considered 
in Eq. (29) reads 



a- (s,r,t) =a% + (2s - 3t)^a^_ . + (4t - 4s)^aL, 



(44) 



with ^Ai^i = a/\T^J as done for the other quantities. As 
one can intuitively expect, the new total uncertainty depends 
on the three free parameters of the model. The resulting models 
are shown in Figure 5. 



4.5 Models M^ and A^s,^ 

The models described in Section 3.5 are derived with a quite 
different approach, which requires an estimate of the mode life- 
time for each star considered in our sample. We note that model 
Mb was also investigated by Huber et al. (2011b), who however 
did not take into account mode lifetimes. 

The natural logarithm of model Mb lead us to 



In 



4(4) 

^bol 



^bol,0 



= -2.5 In 



Z^max,0 



+ 21n( ^1 + 0.51n^ ""^^^ 



+ (2.3-r) In 



'7~osc,0 



(45) 



.^eff,0 , 

Thus, the total uncertainty for the z-th star of the sample is 
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given by 



+ 0.25a? . + (2.3 

'osc,i V 



n2~2 



(46) 



with ^Tosc i = (^Tosc,i/^osc,i^ and is the same for both the models 
here considered, as model A^5,/3 differs only by the additional 
term In (3. The resulting models are shown in Figure 6, with 
similar descriptions as those adopted for Figures 2 and 3. 

4.6 Models Mq and A^6,/3 

Following the same arguments used for the other models, we 
obtain 



l(5) 



In 



^bol,0 



= (2 - 3t) In 
+ (4t - 4) In 



^max,0 



,+0.51n(^^ 

/ \Tosc,Q 



(47) 



+ (4.55 - r - 1.5t) In 



Teff 



,T'eff,0 ^ 

for the model Mq. The new uncertainties to be considered will 
depend on the free parameters r and thus we have 



a1(r, t) = a\. + (2 - 3t)^cr^^ 



, + (4t - 4)'ai,^ 



+ 0.25a? . + (4.55 ■ 



1.5t)^crf^^ 



(48) 



eff,i 



with the same quantities adopted in Eq. (46). An analogous 
discussion to that used for other scaling relations has to be 
applied for model A^e,^- The results of the inference for both 
models are plotted in Figure 7. 



5 MODEL COMPARISON 

As mentioned in Section 4, the term p{A \ M) appearing in 
Eq. (27) (Bayesian evidence) is the one of interest for solving 
the problem of model comparison in the context of Bayesian 
statistics (e. g. see Trotta 2008; Benomar et al. 2009; Handberg 

6 Campante 2011; Gruberbauer et al. 2012). The Bayesian 
evidence is given by integrating the numerator appearing in 
the right-hand side of Eq. (27) over all the possible values of 
the free parameters , ^2 , • • • , ^ A; • Thus we have 



Sm =p{A\M)= [ C{i)i:{i I M)di , 



(49) 



where Q^m represents the parameter space, defined by the in- 
tervals of variation of the free parameters that formalize the hy- 
potheses of the model Al, and having volume given by Eq. (33). 
The Bayesian evidence given by Eq. (49) basically represents 
the integral of the likelihood function averaged by the prior 
distribution. As the prior 7r(^ | M) has to be normalized, 
the evidence depends on the parameter space. Thus, the inter- 
vals [<f j^^^ , <f ™^^] of the free parameters used for computing 
Eq. (49) are those listed in Table 1. 

Since a measure of 8m alone does not carry any meaningful 
information, to solve the problem of model comparison it is 



Table 11. Bayesian Information Criterion (BIC) computed for all 
the models in three cases considered: all targets (second column), SC 
targets (third column), LC targets (fourth column). 



Model 


BIC 


BIC(SC) 


BIC(LC) 


Ml 


7081 


655 


1612 


Mi^p 


5004 


417 


1497 


M2 


8335 


901 


1552 


Ma 


1920 


209 


-393 




M5 


1859 


172 


723 


Mb,^ 


1589 


173 


632 


Me 


1469 


173 


726 


M6,/3 


1468 


137 


635 



useful to take into account the ratios (or odds) of the evidences, 
namely the so-called Bayes factor, which is given as 

p{A I M^) _ 8m^ 



Bij 



p{A I Mj 



(50) 



In case Bij > 1 the model Mi is the favored one, while con- 
versely if Bij < 1 the model Mj ought to be preferred. The 
resulting natural logarithms of the Bayes factor, which are com- 
puted according to Eq. (50) are listed in Tables 8, 9, and 10, 
for the cases of the entire sample, and SC and LC targets, re- 
spectively. Therefore, if In Bij > the model Mi is preferred 
over Mj and vice versa if In Bij < 0. 

It is sometimes useful to consider so-called Information 
Criteria, which may offer a simpler alternative to the Bayesian 
evidence, whose numerical computation in some cases can be 
very time demanding. In particular, we adopted the Bayesian 
Information Criterion (BIC), also known as Schwarz Informa- 
tion Criterion (Schwarz 1978), which follows from a Gaussian 
approximation to the Bayesian evidence in the limit of a large 
sample size, as it can be represented by our sample of stars 
(AT > 1). Thus, the BIC reads 

BIC = -2An,ax + A:lnAr, (51) 

where k is the number of free parameters of the model consid- 
ered (i. e. the dimension of the corresponding parameter space), 
and N the number of data points. Since Amax is known, the BIC 
can be computed straightforwardly. The resulting values of the 
BIC are listed in Table 11 for the cases of the entire sample 
(second column) and of SC and LC targets separately (third 
and fourth columns, respectively). According to the Occam's 
razor principle on which Bayesian model comparison relies, the 
most eligible model is the one that minimizes the BIC. 

As highlighted by the shaded rows and columns, the model 
A^4,/3 is largely the favored one for all the samples considered 
because its evidence is always greater than those of any other 
model investigated in this work. In addition, the BIC confirms 
the result computed through the evidences. 



6 DISCUSSION 

The analysis described in Section 4 and in Section 5 lead us 
to interesting results about the use of the amplitude scaling 
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Figure 2. Marginal PDFs for the free parameters of the models Aii and where expectation values listed in Table 2 (dashed lines) and 

68.3% Bayesian credible regions (light gray for Aii and dark gray for Aii^js) have been marked. 



Table 8. Natural logarithms of the Bayes factor InBij for each pair of models Aii, A4j as derived by means of Eq. (50) 
for the case of the entire sample. 



Model 


Ml 


Ml,/3 


M2 


M4 


M4,/3 


M5 


M5,/3 


M6 


Ml 


















Ml,/3 


1038.2 
















M2 


-626.7 


-1664.9 














Ma 


2582.3 


1544.1 


3209.0 














3690.6 


2652.4 


4317.3 


1108.3 










Mb 


2613.9 


1575.7 


3240.6 


31.6 


-1076.6 








M5,/3 


2748.4 


1710.2 


3375.2 


166.2 


-942.1 


134.5 






M6 


2809.1 


1770.9 


3435.9 


226.8 


-881.4 


195.2 


60.7 




M6,/3 


2810.0 


1771.8 


3436.7 


227.7 1 




196.1 


61.6 


0.9 



relations in asteroseismology. The main aspects of the work 
presented here can be divided into two groups, whose contents 
we discuss in the following. 



6.1 Results from Bayesian parameter estimation 

The results coming from the inference described in Section 4 
show that: 

(i) for models from A^i to A^4,/3, the expectation values of 
exponent r, which we remind was introduced for converting 
radial velocity amplitudes into photometric ones, differ from 
the value r = 2 adopted for non-adiabatic oscillations (Kjeldsen 
k Bedding 1995; Stello et al. 2011; Samadi et al. 2012). On one 
hand, this outcome is even more evident in the case of RGs 
(namely the LC targets), where we found r :^ 10.5 for model 
M2 (Table 4). Although such a high value for r could be partly 
explained by a very tight anti-correlation with the exponent s 
(see Tables 5 and 7), this seems to support the recent findings 
by Samadi et al. (2012) who suggest that for RGs the non- 
adiabatic effects become significantly important in the driving 
mechanism of solar- like oscillations and that, in general, the 
excitation model is underestimating the true amplitudes. On 
the other hand, in the case of SC targets, which are essentially 
dominated by MS stars, the mean r is not as correlated with 
s as for the LC targets (see Table 6), and its estimated values 
are considerably lower (r ^ 3), and much closer to the value 
r = 2 adopted in previous works (Table 3). This again could 
suggest that solar-like oscillations are more adiabatic in early 
stages of stellar evolution. 



Conversely, r becomes negative for the models that take into 
account the granulation power and the lifetime of the modes 
{Mb to Me, 13), which give a different contribution to the pre- 
dicted amplitudes from effective temperature (see Sections 3.5 
and 3.6). The reason our estimates of r are much lower than 
the theoretical value is mainly related to the fact that predicted 
amplitudes based on the scaling relation by Kjeldsen & Bed- 
ding (2011) are considerably greater than the observed ones, 
even for the most adiabatic case corresponding to r = 1.5 (e. g. 
see the discussions by Stello et al. 2011; Huber et al. 2011b). 
This overestimation could be explained by a missing term con- 
taining the mode masses, as proved by recent theoretical cal- 
culations by Samadi et al. (2012) for a sample of RGs observed 
by CoRoT. 

Correlations among the free parameters are very tight in 
many cases and can be partially responsible of decreasing r, 
except for the case of SC targets, where no significant correla- 
tions are found 

(ii) for the models that include the exponent s {Mi to 
Al4,/3), the expectation values derived here are fairly compat- 
ible with previous results (e. g. see Kjeldsen & Bedding 1995; 
Gilliland 2008; Dziembowski & Soszynski 2010; Stello et al. 
2010; Verner et al. 2011; Baudin et al. 2011; Stello et al. 2011; 
Huber et al. 2011b). However, it is worth mentioning that when 
moving from the main sequence to the red giant phase, s de- 
creases, which is apparent for all the models mentioned above. 
This effect can be partially explained by a compensation of the 
exponent r 

(iii) the expectation values of the exponent t found for mod- 
els M4, Me, and Me^/s are not far from the value t = 1.7 =b 0.1 
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Table 9. Same description as for Table 8 but in the case of the sample of SC targets. 



Model 


Ml 




M2 


M4 


M4,/3 


M5 


M5,/3 Me 


Ml 


















118.8 














M2 


-122.4 


-241.2 












M4 


225.2 


106.4 


347.6 












334.7 


215.8 


457.0 


109.5 








M5 


243.7 


124.8 


366.0 


18.5 


-91.0 








242.7 


123.9 


365.1 


17.5 


-91.9 


-0.9 




M6 


244.3 


125.5 


366.7 


19.1 


-90.3 


0.7 


1.6 


M6,/3 


261.4 


142.6 


383.8 


36.2 




17.8 


18.7 17.1 



Table 10. Same description as for Table 8 but in the case of the sample of LC targets. 



Model 


Ml 


Mi,l3 


M2 


M4 


M4,^ 


M5 




Me 


Ml 


















Mi^p 


59.2 
















M2 


30.2 


-29.0 














Ma 


1003.2 


943.9 


973.0 












M4,^ 


1094.2 


1035.0 


1064.0 


91.1 










M5 


445.7 


386.4 


415.5 


-557.5 


-648.6 








M5,/3 


491.5 


432.3 


461.3 


-511.6 


-602.7 


45.9 






Me 


444.9 


385.6 


414.7 


-558.3 


-649.4 


-0.8 


-46.7 




Me,/3 


490.9 


431.6 


460.7 


-512.3 


-603.4 


45.2 


-0.7 


46.0 



derived by Stello et al. (2011) using a sample of cluster RGs, 
although correlation effects with the exponents s and r cannot 
be neglected and are, in some cases, quite pronounced (see Ta- 
ble 5). For model Ma,^ instead, we found t to be very close 
to the value t = 1.32 i 0.02 derived by Huber et al. (2011b), 
who used a very similar sample of stars but with effective tem- 
peratures from KIC. These results are confirmed for both the 
entire sample and the sample of SC targets. The RGs sample 
is instead behaving differently, showing values of t significantly 
lower than the other two samples. This result is certainly dom- 
inated by correlations, which are stronger than those of the 
MS-dominated sample (see below) 

(iv) as a qualitative result we can state that, in general, all 
the expectations of the free parameters involved in the ampli- 
tude scaling relations investigated here are on average rather 
different when comparing the MS-dominated sample to the RGs 
one. The correlations among the free parameters are enhanc- 
ing the differences between dwarfs and giants and are likely to 
be stronger in the sample of more evolved stars because for red 
giants there is a large degree of degeneracy in the stellar funda- 
mental properties as stars with different mass all converge along 
the red giant branch spanning only a small temperature range. 
As a consequence, this could suggest that not a single of the 
amplitude scaling relations discussed can be adopted to both 
MS and RGs simultaneously because the driving and damping 
mechanisms responsible of generating the observed amplitudes 
encounter a substantial change as the stars evolve 

(v) the introduction of a set of models that take into account 
the offset In /3 allowed us to produce new outcomes that provide 
better fits for the entire set of scaling relations adopted, as 
it appears clear by looking at the comparison of different fits 



shown in Figures 3,5, 6, and 7, together with the hint provided 
by our estimates of a weighted rms of the residuals listed in 
Tables 2, 3, and 4. In fact, in almost all cases, the offset ln/3 
differs from zero significantly. In addition, we note that our 
estimates of crj^ns sire fairly consistent with the total relative 
uncertainty adopted in Eq. (29) for models M4, M.4,/3, M.e, and 
Me,^, which are the only ones that have a separate dependence 
upon the mass of the stars (accounted for in the additional free 
parameter t). Models Mi to M2 and models M5 and M^^/s 
have instead a scatter in the residuals that is from 1.5 to 2 times 
larger than the total relative uncertainties. This suggests that 
the average relative uncertainty on the observed amplitudes 
alone is considerably smaller than the intrinsic scatter of the 
residuals given by Eq. (30), for any of the models investigated. 
Thus, we suppose that additional contributions to the total 
uncertainty are missing for models that do not take into account 
a separate dependence upon the mass of the stars. It is however 
important to note that this outcome only holds for the set of 
measurements and uncertainties adopted in the inference. 

6.2 Results from Bayesian model comparison 

From the model comparison applied to the amplitude scaling 
relations here investigated, we can say that: 

(i) models Mi, Mi, 13 and M2 are not performing well at 
predicting the amplitudes along the entire range of stars consid- 
ered. This is clear from the fact that their corresponding Bayes 
factors are the lowest among the models considered (see Ta- 
bles 8, 9, 10), a result that is confirmed by our computation of 
the BIG from Eq. (51), which reaches its highest value for these 
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models (see Table 11 in Section 5). The larger scatter of the 
residuals arising from these models (up to a factor of 2) com- 
pared to the total relative uncertainties derived from Eqs. (39) 
and (42) are supporting our conclusion from the model selec- 
tion described in Section 5 

(ii) models from A44 to Ale,^, according to Bayesian prin- 
ciples of model selection, are certainly more favored than the 
others (Tables 8, 9, 10, and 11). In particular, model A^4,^, 
which includes a separate dependence upon the mass of the 
stars (see Section 3.4), is strongly dominant over all the other 
scaling relations investigated, by at least a factor of ~ exp(880) 
if one takes into account the Bayes factor. Thus, model M.4,/3 
is the best one according to the computation of both the evi- 
dences and the BIC for all the samples used. This suggests that 
the spread observed in the amplitudes is likely to be caused by 
a spread in mass of the stars. This is also confirmed by the 
consistency between crj^ns and the total relative uncertainties 
computed from Eq. (44). 

The models having an evidence weaker than that of A^4,/3, 
but still stronger than that of models from A4i to M2, are the 
models from M.5 to Aie,^, which instead include the effects of 
mode lifetimes and granulation power (see Section 3.5). These 
results are again confirmed for all the three samples considered 
from both the evidences and the BIC 

(iii) all the models that take the proportionality term ^ into 
account, are preferred to their counterparts without this ex- 
tra free parameter, with the only exception of models M.6 and 
Ade,^ for which the model comparison is inconclusive in the 
case of the entire sample of stars (Table 8). These results are 
confirmed from the computation of both the evidences and the 
BIC regardless of the evolutionary state of the stars. As a con- 
sequence, according to the set of measurements and uncertain- 
ties adopted in this work, we can conclude that the Sun is not 
necessarily a good reference star for this Kepler sample, and 
that only models having /3 ^ should be considered. Equiva- 
lently, one could consider models that do not include as a free 
parameter, but having a reference amplitude for scaling given 
by /3Aboi,0 (for models A^2 to M.6,/3) or /3Ax^q (for models 
Ml and Ali,^), with /3 given according to the results shown 
in Tables 2, 3, and 4. Such result could also be supported by 
the effect of the stellar activity (occurring mostly for dwarfs), 
which is responsible of reducing the amplitudes of solar-like os- 
cillations (see the discussion by Chaplin et al. 2011c; Huber et 
al. 2011b). In fact, as noted by Huber et al. (2011b), the Sun 
shows on average a higher activity level than the other stars of 
the sample. 



7 CONCLUSIONS 

A Bayesian approach to test the set of amplitude scaling rela- 
tions discussed in this work allowed us to draw new interesting 
conclusions on the free parameters that describe all the models. 

First, as evident from the results derived in Section 5, our 
analysis strongly recommends the use of Eq. (19) for predicting 
amplitudes of solar-like oscillations, for all the stars spanning 
from MS to RGs. This result, together with the consistency 



found between our estimates of the fit quality, crj^iis? arising 
from the corresponding models and the total relative uncer- 
tainties discussed in Section 4.4, supports the idea that a mass 
difference from star to star is among the main effects that pro- 
duce the observed spread in the oscillation amplitudes (see also 
Huber et al. 2010, 2011b; Stello et al. 2011). 

Second, one should keep in mind that according to the 
Bayesian parameter estimation described in Section 4 the re- 
sults arising from the inference change considerably from MS 
to RGs. This suggests that, in general, particular attention has 
to be paid when using amplitude scaling relations for samples 
containing a large range of stars. This behavior is already ap- 
parent from our plot of the asteroseismic HR diagram (Figure 1, 
bottom panel), which shows that, assuming the logarithmic am- 
plitudes change linearly with the logarithm of the temperatures 
— as done by all the models considered — the two samples have 
on average a quite different slope, which is reflected mainly in 
a different value of the exponent r of the two samples. As a 
consequence, we strongly recommend further investigations of 
this aspect for different sets of measurements and uncertainties 
and also for samples of stars different than the ones used here. 

Moreover, as shown in Tables (5), (6), and (7), the free 
parameters of the models are, in most of the cases, highly cor- 
related. This may suggest to adopt different priors (e. g. non- 
uniform priors) for some of them in future analyses. The corre- 
lation is much more evident for RGs, possibly due to the larger 
range of the fundamental stellar properties in the LC sample 
considered. 

Finally, as it is evident from many of the models investi- 
gated, the Sun is not necessarily a good average star from which 
to scale amplitudes for this sample of Kepler stars, neither for 
the MS-dominated sample nor for the RGs. This outcome could 
be partially explained by the fact that the Sun is placed at the 
edge of our sample of measurements, as visible from Figure 1 
(see also the discussion in Section 4.2), and that its activity 
level is on average higher than that of the other stars (Chaplin 
et al. 2011c; Huber et al. 2011b). Nevertheless, we stress that 
this result also relies on our choice of the set of uncertainties 
adopted for computing the likelihood. In fact, as discussed at 
the end of Section 6.1, the uncertainties used are underesti- 
mated in most of the cases, even when considering errors on 
the variables of the models. This aspect may in fact produce a 
misleading solution to the inference problem and to the model 
selection. Since both the inference and the model comparison 
presented here rely on the assumption that these uncertain- 
ties are themselves error-free, in future work, it will be worth 
investigating how uncertainties affect the results by perform- 
ing uncertainty-independent analyses that do not necessarily 
assume the adoption of a Gaussian likelihood function. 
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Figure 3. Top panels: Predicted amplitudes (filled light-gray circles for LC targets and filled dark-gray squares for SC targets) for model A4i 
plotted against i^max in the three cases considered (all sample in left panel, SC targets only in middle panel, LC targets only in right panel). 
The expectation values of the free parameters reported in Tables 2, 3, 4, have been adopted for plotting the predicted amplitudes. Observed 
amplitudes are shown in background for both SC (open orange squares) and LC (open blue circles) targets, together with 1-a error bars shown 
on both quantities. The Sun's symbol (0) is added for comparison, where Aeso,© = 3.98 ppm (A = 650 nm). The residuals arising from the 
difference between the logarithms of observed and predicted amplitudes, according to Eqs. (38) and (40), are also plotted (same symbols). 
Bottom panels: Same description of the top panels but for the model 




Figure 4. Same description as for Figure 3 but for the model 7^2. 
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